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Abstract 

From the computation of three-point singlet correlators in the two-matrix 
model, we obtain an explicit expression for the macroscopic three-loop am- 
plitudes having boundary lengths (i = 1 ~ 3) in the case of the unitary 
series {p,q) = {m + l,m) coupled to two-dimensional gravity. The sum ap- 
pearing in this expression is found to conform to the structure of the CFT 
fusion rules while the summand factorizes through a product of three modified 
Bessel functions. We briefly discuss a possible generalization of these features 
to macroscopic n-loop amplitudes. 
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The macroscopic n-loop amplitudes describe fundamental processes of c < 1 non- 
critical strings. Their expression is well-known in the case of pure two-dimensional 
gravity realized as the critical point of the one-matrix model it takes a suggestive 
form that the product of factors appears as external legs and the rest are express- 
ible in terms of the sum of the total lengths. These seem to summarize the lore of 
the noncritical string theory in less than one dimension. No comparable expression 
is known for the case of the unitary minimal models coupled to gravity [0]. In 
the previous letter 0, however, we have obtained the two-loop (annulus) amplitude 
having the boundary lengths £i and £2 from the planar solution of the two-matrix 
model for the case of the unitary series; (p, g) = (m + 1, m) 0] (See |l]], 0, p for the 
case of the one-matrix model). Let us first recall this amplitude: ^' Q 

wihj2)c= , 1, E sin-TT K^iMi^)K,_^iMi,) , (1) 

ma^Tl'^ ll + I2 ^ \ m J rn m 

where M = t^/^ and t is the renormalized cosmological constant (The definition of a 
and related objects will be given later in eq. ([I8|) , ([191) .) f\. The product of the modified 
Bessel functions multiplied by the length MiiK^{Mii) is naturally interpreted as 
external leg factors (wave functions). This, together with the absence of a propagator, 
will constitute integral parts of Feynman-like rules of noncritical c < 1 string 0. It is 
also tempting to regard the integer k summed as a label for diagonal primary states 
of underlying CFT, which is now gravitationally dressed. 

In this letter, we will corroborate these features further: we present an explicit 
expression for macroscopic three-loop amplitudes in the case of the unitary minimal 
models coupled to two-dimensional gravity. We obtain this from the computation of 
correlators consisting of the product of three resolvents in the two-matrix model at 
the m-th critical point. We first state our expression for the three-loop w{ii, £2, h)c 
which is manifestly symmetric under the interchange ii ^ ij: 

w{ii,i2,h)c = — . I n 3 3 Hr 

^ This amplitude is for the case in which the spins on the two-boundaries are oriented in the 
same direction. The case of the opposite direction is obtained by replacing £1 + £2 by £1 + {—)™£2 
and by inserting {—)^ in the summand. 

*We have discussed some correspondence of this amplitude with the one from the continuum 
approach (see also 0). 

^ A similar but distinct formula is seen in ||] . 

^Note that the disk-amplitude is also proportinal to (sin^) ^'^^^ Ki^i/rai^f-)- 
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X sin ^TT sin ^vr sin ^vr 

i-^ m m m 

D3(fcl-l,fc2-l,fc3-l) 

X M^eie2i3 K^_^{M£,) K^^^{Me2) K^_,_,{Mis) . (2) 

m m m 

Here we have denoted by V3 

3 3 
1)3(01,02,03) = I (01,02,03) I < 2(m-2), ^Oi = even, 

i=l i=l 

< (e a.) J = 1-3} (3) 

As one may have expected, the product of three factors (sin ^7r)M£ii^'j^(M£j) have 
appeared. 

It is interesting that the selection rules obtained from P3 for the set of integers 
coincide with the fusion rules of the underlying conformal field theory of the unitary 
minimal series (m+ 1, m) In fact, the fusion rules for diagonal primary fields read 
as 

(0n 0ii 0fcfc) 7^ , (4) 

if and only iii+j > k + 1 and two other permutations and i+j + k (= odd) < 2m — 1 
hold. This set of rules is nothing but 2^3 (i — l,j — l,k — 1). 

This coincidence can be understood a little better by studying the small length 
behavior of eq. (H). By letting Mii << 1 independently, we find that eq. (0) becomes 
proportional to 

(—] J2 M^i['''^i\'l'^ti''^l''-^ , (5) 

^ ^ ^ DaCfci-l.te-l.fcs-l) 

where we have used the formula 
vr 



KJx) 



2 sin(i/7r) 



(6) 



h^)-\ 

The M-dependence of the individual terms in eq. (|^) is M~^^ ^ . This is com- 



pared with the computation done in ref. |]T0| (See ref. from the point of view of 
the continuum approach.) from the generalized Kdv flows for the correlation functions 
of the dressed primaries: 

(0, 0, 0,) = iw (|)(^+P+-3)/2^. _ 
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Under u (= t^/™'), the t-dependence of this expression completely agrees with that of 
eq. @. The selection rules for this expression found in ref. is again summarized 
as V3. 

Our formula tells how the higher order operators (gravitational descendants) in 
addition to the dressed primaries are constrained to obey the selection rules of CFT. 
Let us remark here that the c = 1 limit can be studied from our formula ^ . We have 
checked that, by letting one of the three lengths shrink (M£j << 1), the three-loop 
amplitude reduces to the derivative of the two-loop amplitude with respect to the 
cosmological constant t up to a normalization factor. 

The remainder of this letter is devoted to the derivation of eq. (0) from the two- 
matrix model at the (m + 1, m) critical point. We will also briefly discuss what in 
eq. (H) may be generalized to the expression for macroscopic ra-loop amplitudes. 

To begin with let us consider the connected three-point singlet correlator consist- 
ing of arbitrary analytic functions f{A), g{A), and h{A), where the matrix A is a 
dynamical variable. In the bases of orthogonal polynomials, one can derive 

A''~l 00 00 

(Tr/(i) Tr^(i) Tr/^(i)), = E E E [fi^)U [9iA)U [hiA)]u 

i=0 k=N l=N 

N-1 00 7V-1 

- E E E [fi^)h [h{A)]M [9{A)]u , (8) 

i=0 k=N 1=0 

where = {i\f{X)\k). It is convenient to introduce the following 'classical' 

function in the case only the planar limit is of interest, 

A{z,s,A)=J2z'' A{z)k , s = Ai/N , (9) 

k 

where A{i)k = Ai^^^i. The 'classical' function depends on the bare cosmological 
constant A only through s when we take — > 00 limit. It is, therefore, legitimate to 
introduce A{z, s) = lim A{z, s, A) : 

N^oo 

A{z,s,A) = A{z,s) + 0{1/N) . (10) 

In terms of the 'classical' function, we have 

[f{mN)k = ,^<f^J{A{z,s = A)) + 0{l/N) , (11) 

where [/(yl)](z)fc = [/(A)]j_fc j. In the right-hand sind of eq. (H), the leading terms in 
1/A^ of the first term and those of the second term get cancelled. We have to consider 
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the next leading terms. For any integer e, we obtain 



^ Ae 1 / dz dA{z,s) 



df{A{z,A)) 



s=A 



N 2m J z>'+^ ds 
+ (the part independent of e) + 0(l/iV^) . (12) 



The part independent of e comes from the terms 0{1/N) in eq. ([TT]). The second term 
is responsible for the computation in what follows. Using eq. (|12D and considering 
the terms in eq. (||), we obtain 

(Tr/(i) Tr^?(i) Tr/i(i)), 

^ oo oo oo 

~ 7v ^ ^ ^ 

5i=0 52 =0 6=0 

|(52-5i) ^[/(A)](iV),^+,,+i [g{A)]{N)s-S2 [h{A)]{N).s-s.., 
+ {S + 52 + l) [f{A)]{N)s,+s2+i ^j9{A)]{N)s-S2 [h{A)]{N)_s-s,^, 

+ (S- 5,) [/(A)](iV)5,+,,+i [9{A)]{N)s-S2 §-Jh{A)]{N)^s-s,-i } 

+ 0{1/N') . (13) 



Using eq. ([TlD the three-point function eq. (H) in the planar limit can be expressed 
in terms of the 'classical' function as 

(Tr/(i) Tr^7(i) Tr/i(i)), 
A 1 



dzdz'dz" 



N {2niy 7j^[>|z'|>|z"| 

{ (.-.0-(^-.T ^' ^« " 

+ _ _ ^„y f{A(z))g(A{z'))h' {A{z")) , (14) 

where we set A{z) = A{z,s = A). From the above formula (eq. (|14[)), the three-point 
resolvent in the planar limit is expressed as 

N 1 ^ I ^ I 
ir IT IT 



A \ pi- A P2- A P3- A 
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S dzdz'dz" 

(27rz)'^ J\z\>\z'\>\z"\ 

( z dA{z) 



1 



1 



{z - z>Y{z - z"y OA [pi - A{z)]'^ p2 - A{z') ps - A{z") 



+ 



+ 



dA{z') 



{z - z')\z' - z" f px - A{z) dA [p2 - Aiz')]"^ p3 - A{z") 



dA{z") 



{z - z")\z' - z"Y pi - A{z) p2 - A{z') OA [ps - A{z- 



//^l2 



(15) 



where the contour of z encircles that of z' and similarly the contour of z' encircles 
that of z" . We calculate these contour integrals and find 



A ^ Px-A 



Tr 



1 



Tr 



1 



+ 
+ 



P2- A p-i- A'"" 
\ \ dzx d 
o? \ dA dz\ 
dzi dz2 d 



dz\ 



{zi - Z2f{Zx - z-iY dCi 

Z2 dz2 



dz2 dzs 



0(2 dCs 



dCi dA dz2 
dzi dz2 dz3 d 
dCidQi'dAdI:, 



{Zi - Z2y{z2 - Z^y 8(2 



Z3 



dzs 
dzs 



(zi - Z3y{z2 - Zs)^ 



where 



Zi = exp(2?7cosh^j) 
p^~p^ = aCi = 27]"" coshmOi = aM cosh m6'i , r] = 

A - = -(m + 1)7]^"" = -{m + l)aV = -(^ + l)(aM/2)^ 



(16) 

(17) 
(18) 
(19) 



and p*, A,, denote the critical values of p, A respectively. The right hand side of 
eq. (|16D can be written in a compact form as 



N 1 
A \ p^-A 



Tr 



^ Tr ^ 

P2- A P3- A 

1 /aM\-2-V"^ d d d 



2^m(m + l)a^ 



■F{6i, 02, 03 



where 



F{01, 02, 03 



dCi 0(2 0(3 

sinh(m — l)0i 



(20) 



(cosh 6*1 — cosh 6^2) (cosh 6'i — cosh 6^3) 
+ (1^2) + (1^3) . 



sinh 17101 



(21) 
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Making use of the following formula twice, 

1 / sinh(n — k)a sinh(n — k)P \ 

cosh a — cosh P \ sinh na sinh n/? / 

= V sinh(n - j - i + l)a sinh(n -j-k + i)P ^^2) 

~[ ~[ sinh na sinh n/3 



we find that eq. (|2lD is written as a triple sum where the summand factorizes into 
three factors associated with individual loops: 

F{9i, 02, 63) 

^ nv-A: k gin]^(^^ _ f^^Q^ sinh(m - j - i + 1)62 sinh(m - k - j + i)6s 
^ ~{ ~[ sinhm^i sinhm^2 sinhm^s 

(23) 



Here, iy(Ci, (2,(3)0 = f (Tr^ Tr ^ ^r^)^ and £2,^)0 are related by 
the Laplace transform 



00 fOO 



W^(Ci,C2,C3)c = dh di2 di3e-^'''e-^'''e-^''Mh,i2,i3)c (24) 
Jo JO Jo 

= C[w{i^, £2, i3)c] . 

In 0, we have found the following formula for the inverse Laplace image 
^1,(9 sinhfc6'. Mi . kir 

^ = sm—Kj,iMi) . (25) 

Note that Ky{z) is the modified Bessel function. 

Inverse Laplace transform of eq. (^4]) gives the three-loop amplitude in terms of 
each boundary length. We obtain 



1 /oMX"^"^/™ 

M^i,4,4)c = — , 3 3 hr 



m— 1 m—k k 



X y y y sin ^vr sin iri^i^yr sin in^^^n 

fc=l j=l i=l 

X MH^M^ Km-k(Mi,) K ^-,-r+i (Mi^) K ,^^u-j+r (Mi3) .(26) 

mm m 

Let us write the set of restrictions on the triple sum in eq. (p6D as 

JF(A;, j,i) = {(A;, j,i) I 1 < A; < m - 1, 1 < J < m - A;, 1 < i < A;} . (27) 
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By elementary algebras, one can show 

J^ik,j,i) = v^ih-i, k2-i, h-1) 

= V,{k-1, {k + j-z)-l) . (28) 

This completes the derivation of eq. (Q). 

Finally, let us briefly discuss how the features of the three-loop amplitude we 
have found here may be generalized to the n-loop amplitudes. It is for sure that 
the summand contains a product of ra-external leg factors {sin ^^Tc)MiiK_^_kj^{M£i) 

m 

(i = 1 ~ m) but it must be supplemented by a symmetric polynomial, which we need 
to determine. As for the restrictions on the sum, let us first note that V^, which is 
fusion rules, is rephrased as the condition that a triangle with edges (01,02,03) be 
formed whose circumference is even and less than or equal to 2(m — 2). The n-loop 
amplitudes must respect the fusion rules and these triangles will become building 
blocks for the case of ra-loop amplitudes. The ra-gons made out of n — 2 triangles of 
this type 

n 

Vn{ai, ....,a„) = {(ai, ....,a„) | Jl^i = ^ven, 

1=1 

0< |^Ea,j -a„ j = l~n} (29) 

are relevant. Different shapes of n-gons for given (ai,....,a„) and various different 
divisions of each n-gon are important ingredients of the completely and crossing 
symmetric properties of the amplitudes. Also note that the restriction on the sum in 
the case of the two-loop (annulus) amplitudes is written as 'D2{m — ki, ^2) under the 
convention of the factors M^iii2K-^_ki{Mii) K^_k2_{Mi2)- Despite these prominent 
features, complete determination of the amplitudes still requires computation from 
the two-matrix model. 
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